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^ 1 Introduction 

This paper provides the first known exact general solutions of Painleve's sixth equation (PVI) 
and the exact general solutions of the Navier Stokes equations and Prandtl's boundary layer 
equations. In all cases the exact solutions will be given as infinite series expansions. This 
paper is the result of work conducted over a 20 year period on the exact solution of nonlinear 
partial differential equations and nonlinear ordinary differential equations. 

>' 

In around 1991 the author first conducted experimental trial series solution work using Water- 
loo Maple Software and Reduce on the solution of nonlinear differential equations in association 
with G. Smith of the University of Technology Sydney's School of Mathematics, eventually 
leading to publication of work on the Van der Pol equation [17]. Additional, work for the 
author's doctoral dissertation [18] at the University of Wollongong under the supervision of P. 
Broadbridge and T. Marchant led to a further publication on the exact integration of a reduced 
Fisher's equation, a reduced Blasius equation and the Lorenz model[19]. 
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2 Painleve's Sixth Equation, PVI 

According to Sachdev[20, page 423], the problem undertaken by Gambier and Painleve was 
to determine those equations of the form w" = f(w',w,z), where f is rational in w and w', 
whose solutions (as opposed to the equations themselves) are free from movable critical points 
(branch points and essential singularities). Fifty such equations where shown to have solutions 
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having no movable critical points. Of these 50 equations, six in particular required extremely 
elaborate proofs[20, page 424] to show that their solutions had no movable critical points. The 
solutions to the six equations became known as the Painleve transcendents. 

While noting that the Painleve transcendents are, in general, not expressible in the form of 
classical transcendents (solutions of linear ordinary differential equations, exponential, circular 
and elliptic functions) [6, see the definition of classical transcendents at page 330], Ince sug- 
gested the 'possible' existence[6, page 345], of special cases that could be expressible in terms of 
classical transcendents. As noted by Sachdev[20, page 436], Ince's conjecture was confirmed in 
the 1960s and early 1970s by breakthrough work of mathematicians in Belorussia[20, pages 436- 
437], [21, pages 476-477, 888], [10]- [14], which lead to the construction of special solutions in 
terms of Airy, Bessel, Weber-Hermite, Whittaker and hyper geometric functions[20, page 442]. 
Additional work in the important area of construction of special solutions was, as noted by 
Sachdev, also undertaken by Fokas and Ablowitz[20, page 452]. 

Ordinary differential equations whose solutions are of the P type (having no movable criti- 
cal points) have also been seen to have an important relationship to physical properties [20, 
pages 427-429] of waves, in which waves satisfying partial differential equations used to pro- 
duce the ordinary differential equations interact in such a way that their solutions emerge, 
after interaction, substantially unchanged. 

Considerable advances have been made in relation to construction of the exact solutions of 
PI-PVI. Davis produced a very detailed account of the solutions (both analytic and Laurent) 
of PI and PII[5, chapter 8]. Davis further extended his results by use of the method of an- 
alytic continuation [5, chapter 9] and provided methods for determination of the location of 
poles by equating Laurent and analytic expansions in the neighborhood of a singular point[5, 
pages 229-232]. Davis developed the method of analytic continuation both directly (into the 
complex plane) [5, pages 263-266 ] and also by way of the method of "pole vaulting" [5, page 245] 
(in which values of y and y' were found on alternate sides of a pole using a Laurent expansion 
after the location of a pole has been determined [5, page 239-232]. 

Despite an extensive analysis of PI and PII[5, chapters 7-9] using Taylor series and Laurent 
series expansions and also the generation of special solutions for PII-PVI that take the form 
of classical transcendents [10]- [14], no exact general solution of PVI, outside the solutions doc- 
umented in this paper, appear to have been found. Where a solution w = (j)(z, C\,C<i) of 
w" = f(w', w, z) is said to be the general solution[20, page 453] if for appropriate choices of C\ 
and C2, w can be shown to be the solution of an initial value problem in the neighborhood of zq. 



The form of PVI as given by Murphy[15, page 183] has been designated for study. PVI can be 



3 The Exact General Analytic Solution of PVI 



written as 



2x 2 y (x - l) 2 (y - 1) (y - x) y 
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x 2 (x - l) 2 [3y 2 -2y(l + x) + x] p 2 
-2xy (x - 1) (y - 1) [(2x - 1) y + x 2 ] p 

+2 ay 2 (y - l) 2 (y - x f + fix (y - l) 2 (y - x) 2 
+1V 2 {x - 1) (y - x) 2 + bxy 2 (x - 1) (y - l) 2 , 



(1) 
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where in (1) p = y' and q = y" and coefficients ao, «i, a-2 and 03 as appearing in Murphy have 
been replaced by a, (3, 7 and 5 respectively. 

In order to avoid the singularities at and 1 in equation (1) the coordinate system will be 
shifted to —1, by way of the transformation x = x' — 1 after which primes are dropped, and x' 
replaced simply by x. After the transformation x = x' — 1 is effected, equation (1) transforms 
to the following equation: 

y 3 y"[2x 4 - 12x 3 + 26x 2 - 24x + 8] 
+y 2 y"[-2x 5 + 12x 4 - 26x 3 + 24x 2 - 8x] 
+yy"[2x 5 - 14x 4 + 38x 3 - 50x 2 + 32x - 8] 
= y 2 p 2 [3x 4 - 18x 3 + 39x 2 - 36x + 12] 
+yp 2 [-2x 5 + 12x 4 - 26x 3 + 24x 2 - 8x] 
+p 2 [x 5 - 7x 4 + 19x 3 - 25x 2 + Wx - 4] 
+y 3 p[-4x 3 + 18x 2 - 26x + 12] 
+y 2 p[-2x 4 + 14x 3 - 36x 2 + 40x - 16] 

+yp[2x A - 10x 3 + 18x 2 - 14x + 4] (2) 

+2ay 6 

—Aaxy 5 

+y 4 [x 2 (2a + 25) + x(4a + 2/3 + 2 7 - 65) + (-4a - 2/3 - 4 7 + 45)] 
+y 3 [x 2 (-4a - 4/3 - 4 7 - 45) + x{4a + 4/3 + 12 7 + 125) + (-87 - 85)} 
+y 2 [x 3 (2/3 + 2 7 ) + x 2 (2a + 2/3 - 87 + 25) + x(-4a - 8/3 + IO7 - 65) 

+ (2a + 4/3-4 7 + 4<5)] 
+y[-4(3x 3 + 8/3x 2 - 4j3x\ 
+ [2/3x 3 - 6/3x 2 +6/3x - 2/3]. 

The next body of calculations becomes somewhat involved and space does not permit the 
setting out of each element of the calculation. However, part of the first member (the term 
commencing with y 3 y") will be considered in detail, to illustrate the general principles involved. 

Let y = ^2°Z Q aiX % , then 



00 00 

y 2 = a i x% a i x ^ 
i=o j=0 

00 00 

= aia,-**-'". (3) 

i=0 j=0 

The next step in equation (3) is to let i + j = f then 



00 00 

y2 = ^^aiaj/^ixi' . (4) 

i=0 j'=i 

Primes are now dropped in (4) and the order of summation, that is summation by columns is 
used to replace summation by rows. Accordingly equation (4) now becomes 
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oo ] 

y2 = ^2^2 a i a j-i xJ - ( 5 ) 

j=0 i=0 

The practical application of the theory of alternation of summation by rows and columns is 
presented in detail in examples considered in Schwatt[23] together with additional examples 
of operational aspects of series including substantial detail on general expressions for higher 
derivatives of functions and associated Maclaurin series expansions. Additional aspects of the 
theory of infinite series that provide theoretical support for the operational aspects of the the- 
ory of infinite series can be found in Knopp[9] and in Bromwich[2]. Comtet[4] also provides a 
contemporary account of combinatorics and its applications to infinite series. 

The procedure for the development of y 3 will now be given. It will then be obvious that 
the method applied to y 3 can then be systematically developed for higher order differential 
coefficients. From the result for y 2 in equation (5), y 3 can be written as 

oo j oo 

y3 = YY a i a j-i xJ ^2 akxk 

j=0 i=0 k=0 

oo j oo 

= ^2 a i a j-i a k x ^ +k let j + k = k' so k = k! — j 

j=0 i=0 k=0 

oo j oo 

= YjYj 12 a iaj-ia k >-jX k ' 

j=0 i=0 k'=j 
oo j oo 

= a i a j-i a k~j xk on dropping primes 

jr'=0 i=0 k=j 

oo j k 

= ^2 ~^2 12 a i a j-i a k~j xk on changing the order of summation 

k=0 i=0 j=0 

oo k j 

so, y 3 = ^2^2^2a i a j ^ i a k -jX k (6) 

k=0 j=0 i=0 

Next the expression y 3 y" will be fully developed so as to permit equation of coefficients of 
integral powers of x i.e. of x n , 



oo p q oo 

y 3 y" = ^2^2^2ara q - r a p - q x p ^2(s + 2)(s + l)a s+2 x s 

p=0 q=0 r=0 s=0 

oo p q oo 

= ^2 12 12 12^ s + 2 )( s + l)ara q -rap-. q a s+2 x p+s 

p=0 q =0 r=0 s=0 

oo p q oo 

= ^2 12 12 12^ s - P + 2){s - p + l)a r a q - r a p - q a s - p+2 x s 

p=0 q =0 r=0 s=p 

oo s p q 

so, y 3 y" = ^2^2^2^2(s -p + 2)(s -p+l)a r a q - r a p - q a s ^ p+2 x s . (7) 

s=0 P =0 q =0 r=0 
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The full term associated with the differential coefficient y^y" in equation (2) is y 3 y"[2x 4: — 
12x 3 + 26x 2 — 24a; + 8], where the expression 2x 4 in this term can be satisfactorily resolved in 
association with y 3 y" to yield an expression in respect of which coefficients of x % can simply 
be equated. Resolution of the expression y 3 ?/"2x 4 will be considered as representative of the 
calculations involved. From equation (7) we have, after relabeling dummy indices in (7) as 
follows s— > i, p— > j and q— > k, 

oo i j k 

y z y"2x^ = ^^^^(i-i + 2)(i-i + l)a r a fc _ r a i _ fc a i _ j+2 x i 2x 4 (8) 

i=0 j=0 k=0 r=0 
oo i— 4 j k 

= 2^2^2^2^2(i-2-j)(i-3-j)a r a k - r a j - k ai-2-jX l . (9) 

i=4 j=0 k=0 r=0 

Where the transition from (8) to (9) was effected using the substitution i' = i + 4 and then 
substituting in i = i' — 4 and subsequently dropping primes. 

In respect of subsequent equation of coefficients of powers of x l clearly equation (9) is only 
defined for i > 4. This is an important point as some subsidiary equations for coefficients 
are defined only for a fixed value of i and others are defined over a range of values of i, as 
is customarily observed in calculating solutions to linear ordinary differential equations. The 
members found from resolving each member of equation (2) in a manner similar to that under- 
taken above for y 3 y"2x A so as to permit equation of coefficients of x % are listed below. Members 
(1)-(16) below arise from the left hand side of equation (2) and members (17) -(65) arise from 
the right side of equation (2). Upon equating members below, six defining equations need to 
be satisfied, five from i = to i = 4 inclusive and the general equation for i > 5. 
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Members (terms) Arising From Equating Coefficients of x 1 
After Substitution of y = £°^ a 'i x% Equation (2) 



+ 2 Epo ELo £r=o(* - 2 - - 3 - i)a r a fc _ r a i _ fe a i _ 2 -j 

-12 Y!j~=Q ££=0 £r=o(* - 1 - J')^ - 2 - j)a r a k - r aj- k ai-].-j 

+26 ELo £r=o(* - - 1 - j)ara k -raj-k(H-j 

- 24 Ej=0 ELo £r= (* + 1 - - fiarak-raj-kdi+i-j 

+ 8 £J=o ELo £r=o(* - J + 2 )(* - J + l)a r .a jfc _ r a j _ fe a i _ J+2 
~ 2 Epo £Lo(* - 3 - - 4 - fidkaj-kai-z-j 
+ 12 Epo £Lo(* - 2 - - 3 - j)akaj-kai-2-j 
-26 Epo £fc=o(* - 1 - - 2 - j)a>ka>j-kai-i-j 
+ 24 £Lo Efe=o(* ~ - 1 - j)akaj- k ai-j 



(* + 1 - - j)akaj-kai+i- 



+ 2 £}=o(* -3-j)(i-4- j)*j<H-z-j 
-14Ej=o(* - 2 - - 3 - i)a i a i _ 2 -j 
+38Epo(» - 1 - - 2 - jjaj-oi-i-j 
-50E}=o(* ~ - 1 - j) a i a i-j 
+ 32 E}=o(* + 1 - - j)o-jO-i+l-j 

~ 8 EJ-=o(* + 2 ~ •?)(* + 1 ~ j) a 3 a i+2-j 

+ 3 E}=o ELo Er=o~ J ( r + 3 - i - r)a fe a j _ fe a r . + ia i _ 3 _ J _,, 

-18 Epo ELo Er=o~ J ( r + 1)(« - 2 - J - r)a fc o i _ fc a r+ iai_ 2 -j- 
+39E}=o ELo Er=o~ J ( r + 1)(« - 1 - J - r)a k aj-kar+iai-i-j- 
-36E}=o ELo Er=o~ J ( r + !)(*- J - r)a fc a j _ fc a r+ iai_ i _ r 
+ 12 E}=o ELo Er=o( r + + 1 - i - J")a fc a J -_ fc a r+ ia i+ i_j_ r . 
~ 2 £po Ei=o( /c + 1)0' - fc + l)afe+iaj-fc+iai-5-i 

+i 2 E}=oELo(^ + i)(j-^ + iH+ 

-26 Epo Efe=o( A; + - fc + lH+ia.j-fc+iai-3-j 
+ 24 £po Ei=o( fe + - + l)a fe+ ia J -_ fc+ ia i _2-j 
- 8 E}=oEi=o( /c + ~ fc + l) a fc+i a j-fc+i a i-i-j 
+ E}=oO' + - 4 - j')a j+ iai_ 4 -j 
- 7 Ej=oO' + ~ 3 ~ i) a i+i a i-3-j 
+19E$=o(i + " 2 " i)a j +ia i _ 2 - j 
- 2 5EpoO + !)(* - 1 - ./)"/ • i "' : j 
+16 E}=o(i + !)(*- j)aj+iai-j 
- 4 E}=o(i + !)(*- i + l)aj+iai_j + i 
~ 4 Epo ELo £r=o(* ~ 2 ~ j)a r a k - r aj-kai-2-j 
+18Epo ELo Er=o(* - 1 - j)a r ak- r aj-kai-i-j 
-26 Epo ELo Er=o(* - J>rOfc-raj-fcai-i 
+ 12 E}=o ELo Er=o(* - J + l)a r a k - r a j - k a i - j+1 
~ 2 E}=o £k=o(* - 3 - i)afcai-fcOi-3-j 
+ 14 Epo ELo^ ~ 2 ~ j) a fc a j-fc a *-2-j 
-36EpoELo(^ - 1 - J>fcOj-fcai-i-j 

+ 40 Epo Ei= ( i - i)«fc a i 

-!6E}=o £fc=o(* - + l)afcOj-fcOi-j+i 



■r 1 — 



i > 4 


[1] 


i > 3 


m 


i > 2 


[3] 


i > 1 


t 4 ] 


i > 


[5] 


i > 5 


[6] 


i > 4 


[7] 


i > 3 


[8] 


i > 2 


[9] 


i > 1 


[10] 


i > 5 


[11] 


i > 4 


[12] 


i > 3 


[13] 


i > 2 


[14] 


i > 1 


[15] 


i > 


[16] 


% > 4 


[17] 


i > 3 


[18] 


% > 2 


[19] 


i > 1 


[20] 


i > 


[21] 


i > 5 


[22] 


i > 4 


[23] 


i > 3 


[24] 


% > 2 


[25] 


i > 1 


[26] 


i > 5 


[27] 


i > 4 


[28] 


i > 3 


[29] 


i > 2 


[30] 


i > 1 


[31] 


i > 


[32] 


i > 3 


[33] 


% > 2 


[34] 


i > 1 


[35] 


i > 


[36] 


i > 4 


[37] 


i > 3 


[38] 


% > 2 


[39] 


i > 1 


[40] 


i > 


[41] 
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+2 y^-_ n (i — 3 — j)aj(ii-3-j 


i > 4 


[42] 


— 10 y^Ln(« — 2 - j)ajai-2-j 


i > 3 


[43] 


+ 18 E}=o( i ~ 1 ~ j) a j a i-l-j 


i > 2 


[44] 


-14 YZ-i \(i -j)ajai-j 


i > 1 


[45] 


+ 4 Ei-o(* + 1 ~ j) a j a i+l-j 


i > 


[46] 


+2a Yjk-ft Er- (1 EI-O Ef-(1 ^tO-s-tO-r-sO-k-r^j-k^i-j 


i > 


[47] 


—4a Ej=0 Efc=0 Er=0 Es=0 a s^r-sO'k-r a j~k a i-l-j 

+(2a + 25) EPo Efc-n Er-o a rak- r aj~kai-2-j 


i > 1 


[48] 


i > 2 


[49] 


+(4a + 2/3 + 27 - 65) EPo E£-n Er-o a r ak- r aj-kai-i-j 


« > 1 


[50] 


+(-4a - 2/3 - 47 + 45) E1-0 El- n Er-n a r a k - r a,j- k ai-j 


i > 


[51] 


+(-4a - 4/3 - 4 7 - 45) Et-o El-n o fc a 7 -_ fc ai_ 2 -j 
+(4a + 4/3 + 127 + 125) Et-n Ei n a fc a,_ fe ai_i_j 


i > 2 


[52] 


i > 1 


[53] 


+ ( — 87 — 85) y^!- n n akCtj-bClj-j 


i > 


[54] 

L J 


+(28 + 2t) a.- a,- ^ „■ 


i > 3 


[551 


T^a + — 07 -f- AO) Z^i j=0 a ] a i-2-j 






+(-4a - 8/3 + IO7 - 65) }^j=o a j a i~i-j 


Z > 1 


[57J 


+(2a + 4/3 - 4 7 + 45) ^ i=0 


i > 


[58] 


-4/3aj_ 3 


i > 3 


[59] 


+8/3ai_ 2 


i > 2 


[60] 


-4/3aj_i 


i > 1 


[61] 


+2/3 


i = 3 


[62] 


-6/3 


i = 2 


[63] 


+6/3 


i = 1 


[64] 


-2/3 


i = 


[65] 



Next a number of simple observations can be made about members 1-65, hereinafter referred 
to as M{ for i = 1 . . . 65. It is to be noted that the highest possible subscripts can only 
be obtained from members 5 and 16 by setting j = 0. This fact permits the separation of 
member 5, M5, and member 16, Miq, in set (10) above, into two separate parts so as to put 
the defining equation for the coefficients arising from substitution of y = a * x * ^ n ^° (^) 

into a form that can be used not merely to define the structure of the solution to equation (2) 
but also to enable the solution to be constructed iteratively. This is achieved in detail as follows. 

Members 5, (M5) can be expanded as follows 

8 E}=o Ei=o Er=o arak-raj-kai-j + 2{i -j + 2)(i-j + 1) 

8ala i+2 (i + 2)(i + 1) + 8^- = i EiLo Ylr=o a r a k - r aj- k ai- j+2 (i - j + 2)(i - j + 1) 
(i > 0) (i > 1) 

Similarly, member 16, Miq, appearing in set (10) can be expanded as follows 



-8 E}=o ajO-i+2-j(i + 2 - + 1 - j) 

-8a ai +2 (i + 2)(i + 1) - 8 E} =1 a 
(i > 0) (i > 1) 



(12) 

8a a i+2 (i + 2)(t + 1) - 8 E}=i aiai+2-i(« + 2 - + 1 - i) 
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Members of set (10) can now be summed and solved for cij + 2 after member 5, M5, and member 
16, Mi6, of equation (10) have been expanded in accordance with equations (11) and (12) to 
yield the following defining equation for the coefficients of the solution to equation (2) 



- V 16 M* + V 65 M* 

a *+ 2 = trlov-lur? A for *- 0; ao ^ 1 ' (13) 

8a (t + 2)(i + - 1) 

It is to be noted that the members appearing in equation (13) now have a star over them to 
emphasize that members 5 and 16 (M15), (M±q) have been modified to the extent that sum- 
mation now only occurs from j = 1 onwards. 

To demonstrate the validity of equation (13) parameters and values of ao and a\ have been set 
as follows a = /3 = 7 = (5 = l, ao = 2 and a\ = 1. After, evaluation using equation (13) the 
following values were found. 

a = 2, 01 = 1, a 2 = 5/48, a 3 = 311/864, a 4 = 18725/20736, a 5 = 48313/34560, a 6 = 
17430769/8957952, a 7 = 3838061/1451520, a 8 = 65037559477/18059231232, 
a 9 = 95777442903929/19503969730560. 

These values have been confirmed by using Waterloo Maple Software's dsolve routine. 

Laurent expansions in the neighborhood of a singularity have also been conducted by the au- 
thor. However, as the the construction of these expansions involve the use of substantially no 
new principles beyond those recited above no further elaboration on this point will be under- 
taken. 

The next step will be to show how the methods used above to solve ordinary differential equa- 
tions exactly can be generalized to solve partial differential equations exactly. The general 
system of the Navier Stokes equations will be solved exactly and the methodology derived 
above will then be applied to solve Prandt's boundary layer equations. 



4 The Navier Stokes Equations and Boundary Layer Flow 

Ludwig Prandtl said of the Navier Stokes equations, "The general differential equations of 
viscous fluids lead to mathematical difficulties which may be unconquerable for a long time to 
come." ... ".[16, page 104] 

This paper presents what is apparently the first set of general and exact solutions of the Navier 
Stokes equations as well as presenting what appears to be the first set of general and exact 
solutions to the associated boundary layer equations of Prandtl. 

At large Reynold's number the flow around a streamlined body, with no flow separation, can 
be divided into a boundary layer flow confined to a region of thickness 5 around the body (in 
which viscous effects are important) and a flow field far from the body which closely resembles 
non viscous flow[8, page 160] (where the Reynolds number is given by Re = pUL/p = UL/u 
for characteristic velocity U, length L and density p and where 5 <C L, coefficient of viscosity 
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[x, and kinematic viscosity v). 



Under extreme conditions the velocity profile around the body can become extremely dis- 
torted until separation and associated eddying and turbulence occurs[8, pages 160 and 168]. 
In principle separation is characterized [22, page 113] by {du/dy) y= Q = 0. However, in practise 
separation is far more complex and turbulence can even occur well before separation, leading 
paradoxically to the flow adhering[8, page 161] to the body wall even in the presence of an 
adverse pressure gradient. 

Once the boundary layer equations have been solved "some" estimate of the location of the 
point of separation can be determined. However, it is to be noted that in practise, as recited 
above, despite 'theoretical' calculations as to the point of separation, other factors come into 
play, including the fact that the presence of any turbulence in the flow could actually delay 
flow separation. Accordingly, results from theoretical calculations must used in association 
with experimental results to select appropriate modifications to body structures (typically air- 
craft wing and ship hull structures) that minimize drag and in the case of an aircraft wing, 
increase lift. 

Typically, the modifications involved in boundary layer flow stabilization include the introduc- 
tion of slats, slots and slits so as to delay the onset of flow separation [2 2, page 268]. A detailed 
account of which is given in Schlichting under the heading "Boundary Layer Control" [22, 
chapter XIII]. By way of comment on flow stabilization techniques, it is to be noted that even 
very slight adjustments to body structure by way of the introduction of slits, slats and other 
structural changes designed to maintain stable laminar flow, can yield extremely significant 
increases [22, page 275] in lift and also significantly reduce drag. 

The Navier Stokes equations and the associated boundary layer equations will now be solved. 



5 The Exact General Solution of the Navier Stokes Equations 

The Navier Stokes Equations are written as follows (on the assumption that there is no free 
surface and that pressure is taken relative to static pressure were no flow to exist). The ab- 
sence of the body force term also implicitly assumes that no body force other than gravitational 
forces exist, an assumption that is reasonable for most flow structures likely to be encountered 
outside highly specialized flows that arise in areas including magneto-hydrodynamics. 



Equation (14) is to be solved in conjunction with the continuity equation so as to provide four 
equations in four unknowns[7, page 144] u, v, w and P. 

The continuity equation is given by 




P 
1 

P 



(14) 
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No new principles of solution methodology are introduced beyond those already considered in 
the solution of PVI. However, for the sake of completeness, the general principle in relation to 
resolution of nonlinearity for purposes of equating coefficients will again be illustrated in some 
detail. 

For u = ^i=o Sj=o Sfc=o Sz=o ^-ij,k,ix l y^ z k t l the expression udu/dx will be resolved so as to 
permit equating of coefficients. 



r-. r-. OO OO OO OO 

du d 



dx dx 

i=0 j=0 k=0 1=0 



OO OO OO OO 



EEEE^,^~v^ let i = i'+i 

i=0 j=0 k=0 1=0 



du 
dx 



OO OO OO OO 



E E E E^ ^-)Ai+i,j,k,l% l y J 'z h t l after dropping primes (16) 

i=0 j=0 k=0 1=0 



In calculating terms including du/dx above, terms including i = — 1 that result in a value of 
zero associated with i + that arose during use of the transformation i = i' + 1, see (16), have 
obviously not been included. Using (16) we obtain, 



du 
dx 



oooooooo oooooooo 



E E E E A P , q ,r, s x^z r f EEEE(* + m+^wzH 1 

p=0 q=0 r=0 s=0 i=0 j=0 k=0 1=0 



oooooooooooooooo 



EEEEEEEE(* + i)^,^m.,M* i+ v +< ^ + ^ +s 

p=0 q=0 r=0 s=0 i=0 j=0 k=0 1=0 



oooooooooooooooo 



du 
dx 



EEEEEEEE(*-^+ —p+l,j—q,k—r,l—s x ]P Z t 

p=0 q=0 r=0 s=0 i=p j=q k=r l=s 

oo oo oo oo i j k I 

EEEEEEEE( ,; -p- i)^ l , iria ^_ p+ i J _, )fc ^ ) ,_^yz¥ (17) 

i=0 j=0 k=0 1=0 p=0 q=0 r=0 s=0 



Where as before in the construction of equation (17) the transformations i = i' — p, j = j' — q, 
k = k' — r and I = 1' — s and the method of changing the direction of summation from rows to 
columns were used. 

No new principles other than those recited above have been used in the resolution of other 
elements of the Navier Stokes equations. Accordingly, setting 

oo oo oo oo 

^EEEE^.^v'A ! 

i=0 j=0 k=0 1=0 

oo oo oo oo 

^ = EEEEw^ 

i=0 j=0 k=0 1=0 
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and 



oo oo oo oo 



i=0 j=0 k=0 1=0 



then the coefficients of the analytic solution of the Navier Stokes equations (14) are defined 
by the following equations (after resolving each nonlinear member as occurred in relation to 
udu/dx so as to permit equating of coefficients of x l y J ' z k t l ) 



^2p=0 ^2q=0 Sr=0 Ss=0 ~~ P + ^)^P,q,r,s^-i-p+l,j-q,k-r,l-s 

Q ~l~ ^)-^p,q,r,s-^-i— p,j— q+l,k— r,l— s 
+ (k - r + l)C Pi9iriS A-p,j- 9 ,fc-r+l,«-s) + {I + l)Ai,j,k,l+l 

+ l)A+ij,fc,i + "((* + 2)(i + l)A+2j,ifc,/ + (J + 2)0' + 1)A, j+ 2,m 
+(k + 2)(k + l)A id , k+2 ^ 



J2p=0 Y^q=Q Ylr=0 Ss=0 ((* ~~ P + l)-<4p,g,r,s A-p+l,.j-g,fc-r,Z-s 
+ (j — <Z + l)Bp,q,r,sBi—p,j—q+l,k—r,l—s 

+ (k - r + l)Cp )?)r)S Sj_pj_q ife _ J . +li ;_ s ^ + (/ + l)i?i,j,fc,Z+l 

ft? + i)Aj+i,fc,i + + m + i)A+2j,fc,i + (i + 2)(j + i)Aj+2 ) fc 1 j 

+(fc + 2)(fc + l)A Jifc+2)J ) 



(18) 



(19) 



Sp=0 Sg=0 Sr=0 Ss=0 P + ^)Ap,q,r,sCi-p+l,j-q,k-r,l-s 

+ (j — 9 + l)-^p,g,r,sA— p,j— g+l,fc— r,i— s 

+ - r + l)Cp j g jrjS Ci_p j _ 7 -_ gj fe_ r . + i j j_ s ^ + (/ + l)Cjj i A; l / + l 

(zUJ 

+ lJAj.fc+i.J + f ((» + 2)(i + l)C i+2j - ) fc ) j + (j + 2)(j + l)C ij - +2) fc,i 
+(fc + 2)(fc + l)C ij - fc+2)J ). 
The defining equation for the coefficients arising from the continuity equation is 

= (i + l)A i+lij>ki i + (j + l)B itj+1A i + (k + \)C id>k+lj i (21) 

6 The Exact General Solution of the Prandtl's 
Boundary Layer Equations 

A highly simplified form of the general Navier Stokes equations (14) , known as Prandtl's 
boundary layer equations (including the continuity equation) is given as follows [22, page 109] 
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Together with boundary conditions = u = vaty = and u = U(x, t) at y = oo. 



For 



and 



oo oo oo 

^EEE^v^ (24) 

i=0 j=l k=0 



oo oo oo 

v = EEE B ^ (25) 

i=0 j=l k=0 



a simplified equivalent version of the general defining equations for the coefficients of the so- 
lution of the general form of the Navier Stokes equations arises. It being noted that for u and 
v, the starting value for the index j is in fact 1 and not 0, as a result of the fact that u = at 
y = as given above. (For y = we obtain = -B^o.fc = ^i,o,fc as au other terms vanish in u 
and v .) 

Further, as noted by Schlichting as pressure is impressed[22, page 110] upon the boundary 
layer by the external flow, the pressure gradient is fully defined and given by the following 
expression relating to distal flow U (x, t) being 

dU_ JJ dU_ ldP_ 

dt dx p dx 

After substitution of u and v as defined by equations (24) and (25) into Prandlt's boundary 
layer equations (22) and (23) and then using (26) and applying the same methods as used to 
resolve the nonlinear terms in PVI and in the general form of the Navier Stokes equation, the 
following defining equation is obtained 

Er= j2T=ij2k=oji,j,^(k + i)xV'«* 

+ J2i=0 £j=2 Sfc=0 Sp=0 £q=l J2r=0 \,q,r-^i-p+l,j-q,k-r{^ ~ P + l)x l yH k 



+ ZT=o ££i ET=o E;=o T.U 52=0 B^rA^.^rU - q + l)xV't* 

Er=oEr= ^ + i(Hi)rf 

+ E£o Elto E;=o E k q=0 U p , q U w ^ q {i -p + l)xH k 
v ET=o Er =0 Ello Ai, j+ 2, k (j + 2)(j + i)*V'f* 



where tf(s,t) = 5Xo 



(27) 
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Equation (27) can be further resolved by noting that only under the condition j > 2 does 
a general defining equation for the coefficients result. Accordingly two special subsidiary 
equations result for the special cases of j = and j = 1 being 

i k 

= U itk+1 (k + 1) + E U P , q Ui- p+ltk - q {i -p + 1) + 2uA iXk (28) 

p=0 q=0 

for j = and i, k > 
and 



+ i ) + e e b p 

,l,r-^-i—p,l,k—r 

p=0 r=0 

which can also be written as 



Ai, 3 ,k = ( fc + + E E 5p,l,r A-p,l,fc-r] (29) 

p=0r=0 

for j = 1 and i,k >0. 

Equations (28) and (29) define -Aj 2,fe an d 3 ^ respectively. Similarly, the defining equation 
for A i: j + 2,k for j > 2 and i,k > can be extracted, by simple equation of coefficients, from 
(27) to yield 

Aj+2,k = f . n( , m (Xj,fc+l(& + 1) + E E E(* " P + l)Ap,g,rAi_p + i J _ gifc _ r 

' p=0 9=1 r=0 

i j k 

+EEE 

^p,q,r A i-p,j-q+l,k-r (j -9 + 1) (30) 

p=0 <jr=l r=0 

for j > 2, i, k > 0. 

Equations (28), (29) and (30) can be further simplified using relations derived from the conti- 
nuity equation as follows. 

By substituting u = ££=1 Efe==o A i,j,kX % y 3 t k and w = £)~ £~ x YX=o B i,j,kX l y 3 t k into 

the continuity equation (23) the following equation is obtained 

OO OO OO OO OO OO 

= E E E A^uA* + W + E E E ^ + i W* fc . (31) 

i=0 j=l fc=0 i=0 j=0 k=0 

Upon setting j = in (31) we obtain- 
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= B iAjk (32) 
Further, the general defining equation from equation (31) for i, k > 0, j > 1 is 

= Ai +1Jik (i + 1) + B iJ+hk (j + 1) 
which can be written in a slightly more useful form, for computational purposes, as 

Bi,j+i,k = — '. _j_ ^ Ai+i,j,k (33) 
From equation (32) the defining equation for A ij3 ^, (29) becomes simply- 

Ai A k = -^[Ai tl ,k+i(k + l)] (34) 

for i, k > 0. 

From equation (33) we obtain 

ft _ P +1 A 



which when substituted into equation (30) yields 

A ^+2, k = v(i + l )(i + 2) (A,j,k + i(k + 1) 



Ki + i)(j + 2) 

i j—1 k 

+ ~ P+ ^Ap^rAi-p+lj-q^k-r (35) 

p=0 q=l r=0 

1 3 * (p+l)(j_- g +l) 
<7 



Vyy (p+l)Cj-g + l) , A \ 

2_^i „ ^p+l,g-l,r^i-pj-g+l,fc-r^ 



=0 g=2 r=0 



for i,k > and j > 2. 



Not only are the coefficients in the expansion of u and v now fully defined but the coefficients 
can also be systematically constructed from the above results. 

The above results can be summarized as follows 

Ai,2,k f° r h k > is fully defined by equation (28). 

^4i,3,fc for i, k > is fully defined by equation (34). 

Aij+2,k f° r h k > j > 2 is fully defined by equation (35). 

From the above it is seen that the only remaining unknown is -Aj i However, Ai i k is fully 
and completely defined by way of matching u against the external flow. In principle the match 
against the external flow occurs in the limit as y — > oo, of course in practise the match occurs 
over a dimension, that is small in comparison with the reference length L. Accordingly, the 
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solution u can be matched against the external flow at a finite distance from the body surface 
so as to define A^^. 

The above process will fully define u. It is also noted that v is also fully defined since J3j j + i ^ = 
— ((i + l)/(j + l))Ai + i : j t k for i,k > 0, j > 1 from equation (33). Accordingly, from equation 
(33) we can write 

Bi,j,k = j—A+i,j-i,k (36) 

for j > 2 which together with equations (28)-(35) fully define A,j,fe (where = -B^o.fc = Bi,i,k 
fori,fc>0). 

A few further simple adjustments will be made so as to put u and v in a form that is more 
suitable for calculation. 

Equation (35) can be written, by letting j = j' — 2 and subsequently dropping primes as 



Ai ^ k = v (j _ (^J-2,fc+l(fe + !) + mi Zfo ~ P + 1 )A»,g,r^«-p+lj-2- g ,fc- 

p=0 g=l r=0 



i j-2 fc 



2-^i 2-^i q /i p+l,q-l,r^-i-p,j-l-q,k-r J {<><) 



(p + l)(j-l-q) 

p=0 q=2 r=0 

for i, k > and j > 4. 



Further, so as to use equations (28), (34) and (37) to best effect for purposes of computation, 
u will be written as 



oo oo oo 
i=0 j=l k=0 



= YjYj A i,hkx l yt k A i,2,kx i y 2 t k + A i,3,kxVt h 

i=0 k=0 i=0 k=0 i=0 k=0 

oo oo oo 

i=0 j=4 k=0 

oo oo oo oo oo 

= 22 [ A ^y + A lXk y 2 + A iAk y^ xH k + A i,J,kX i y j t k 
i=0 fc=0 " i=0 j=4 k=0 

Equation (37) and a similar equation for v derived by using equation (36) can be used to 
construct the solutions to Prandtl's boundary layer equations. These methods have been used 
to generate solutions which have then been substituted back into the boundary layer equations 
by using Waterloo Maple Software to verify the validity of the expansions. It has then been 
confirmed, that the solutions constructed to increasing orders of magnitude, lead to remainders 
(residuals) that systematically vanish as the order of expansion is increased. 
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7 Conclusions 



It is seen that the above methods can be used to construct exact analytic solutions to many 
differential equations (ordinary and partial) that are likely to be encountered in practise. Nev- 
ertheless, many challenges remain for further research. The use of exact analytic expansions 
in numerical analytic continuation as outlined in Davis. Additional challenges of a non-trivial 
nature also include the desirability of constructing special solutions to differential equations in 
terms of classical transcendents (as noted above in relation to the work conducted by mathe- 
maticians in Belorussia (Lukashevich and Yablonski) and also by other contributors including 
Fokas and Ablowitz on Painleve transcendents). 

In closing, it has to be acknowledged that it is completely impossible in a paper of this size to 
do justice or pay tribute to the many important works worthy of further consultation in the 
areas of fluid dynamics and special functions. However, two works, by two leading authorities 
in their fields, provide a starting point, at the very least, for purposes of further investigation. 
In the area of fluid dynamics an excellent classical work that provides a thorough grounding in 
the principles of fluid dynamics is G. K. Batchelor's An Introduction to Fluid Dynamics[l] and 
in the area of the Painleve transcendents, P. A. Clarkson's review of the Painleve transcendents 
in the NIST Handbook of Mathematical Functions^} also provides an excellent starting point 
for further research. 
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